Electronic properties of closed cage nanometer-size spherical graphitic particles 
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• We investigate the localization of charged particles by the image potential of spherical shells, such 

as fullerene buckyballs. These spherical image states exist within surface potentials formed by the 
competition between the attractive image potential and the repulsive centripetal force arising from 
■ the angular motion. The image potential has a power law rather than a logarithmic behavior for 

a nanotube, leading to fundamental differences in the forms for the effective potential for the two 
geometries. The sphere has localized stable states close to its surface. At low temperatures, this 
results in long lifetimes for the image states. We predict the possibility of creating image states 
with binding energies of a few meV around met allic/non- metallic spherical shells by photoionization. 
Applications and related phenomena are discussed. 

PACS numbers: 78.40.Ha, 34.80.Lx, 73.20.Mf,32.80.Fb 
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The experimental and theoretical study of carbon is currently one of the most prevailing research areas in condensed 
matter physics. Forms of carbon include several allotropes such as graphene, graphite as well as the fullerenes, which 
cover any molecule composed entirely of carbon, in the form of a hollow sphere, ellipsoid or tube. Like graphite, 
Q ■ fullerenes are composed of stacked graphene sheets of linked hexagonal rings. For these, the carbon atoms form 
(_J ' strong covalent bonds through hybridized sp 2 atomic orbitals between three nearest neighbors in a planar or nearly 
planar configuration. 

Mass spectrometry experiments showed strong peaks corresponding to molecules with the exact mass of sixty carbon 
atoms and other carbon clusters such as C70, C76, and up to C94 [l|,|2|. Spherical fullerenes, well known as "buckyballs" 
(Ceo), were prepared in 1985 by Kroto, et al. [3| The structure was also identified about five years earlier by Iijima, 
[J] from an electron microscope image, where it formed the core of a multi-shell fullerene or "bucky onion." Since 
then, fullerenes have been found to exist naturally @. More recently, fullerenes have been detected in outer space 
Q. As a matter of fact, the discovery of fullerenes greatly expanded the number of known carbon allotropes, which 
until recently were limited to graphite, diamond, and amorphous carbon such as soot and charcoal. Both buckyballs 
and carbon nanotubes, also referred to as buckytubes have been the focus of intense investigation, for their unique 
chemistry as well as their technological applications in materials science, electronics, and nanotechnology M . 
' Recently, the image states of metallic carbon nanotubes Band double- wall non- metallic nanotubes [J 3 were 
investigated Experimental work [ill ] includes photoionization |12| and time-resolved photoimaging of image-potential 
states in carbon nanotubes [l3| . There has been general interest in these structures because of electronic control 
on the nanoscale using image states. This has led to wide-ranging potential applications including field ionization 
of cold atoms near carbon nanotubes, [l5[ and chemisorption of fluorine atoms on the surface of carbon nanotubes 
[l6l |. Here, we calculate the nature of the image-potential states in a spherical electron gas (SEG) confined to the 
surface of a buckyball in a similar fashion as in the case of a nanotube. For the semi-infinite metal/vacuum interface, 
[TtJ related image-potential states have been given a considerable amount of theoretical attention over the years. 
Additionally, these states have been observed for pyrolytic graphite 18] and metal supported graphene 19]. Silkin, et 
al. (20| further highlighted the importance of image states in graphene 21!] by concluding that the inter-layer state in 
graphite is formed by the hybridization of the lowest image-potential state in graphene in a similar way as occurs in 
bilayer graphene [22, 23|. The significance of the role played by the image potential has led to the observation that for 
planar layered materials, strongly dispersive inter-layer states are in common. However, the eigenstates for a spherical 
shell are non dispersive and so too are the collective plasma modes, f2~j - l2"8j leading to interesting properties for the 
image potential. Fullerene is an unusual reactant in many organic reactions such as the Bingel reaction that allows 
the attachment of extensions to the fullerene [29] . For this reaction, the shorter bonds located at the intersection of 
two hexagons (6-6 bonds) are the preferred double bonds on the fullerene surface. The driving force is the release of 
steric strain. The surface-related properties of these structures are of course influenced by its image potential, making 
it one reason why we are interested in this potential. 

The fullerenes are modeled in analogy with the collective excitations in a planar semiconductor two-dimensional- 
electron-gas. We now consider a spherical shell of radii R whose center is at the origin. The background dielectric 
constant is ei for < r < R, and €2 for r > R. An electron gas is confined to the surface of the sphere. If a charge 
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FIG. 1: (Color online) Schematic illustration of a charged particle captured by the image potential and orbiting around a 
buckyball. The radius of the orbit is determined by the dielectric constant within and surrounding the shell as well as the 
angular momentum quantum number of the captured particle. For semiconducting shells, the localization is strong and the 
radius of the stable orbit can be a few nanometers. The localization is weak for metallic shells. 

Q is located at (ro,0o,4>o) in spherical coordinates, then for tq > R, the total electrostatic potential is given by 
$tot = ^cxt + ^md, where $ e xt is the external potential due to the point particle and $i n d is the induced potential. 
The external potential may be expanded in the form 



$ cxt (r, 9, 0) = AvkQ ]T 1 %y i M ( fi o) Ylm (fi) , 

LM r > 
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where k = (47reo) 1 with eo the permittivity of free space. Also, Ylm (O) is a spherical harmonic and 51 is a solid 
angle. For r < R we express the total potential in the form 



*£> (r, 9, 0) = AvkQ 77f^-7 A L r L Y LM (O) 

LM ^ + 1 
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whereas for r > R we express the induced potential as 



(r, 0, 0) = AvkQ J2 7^-lB L r-^Y LM (fi) 

LM 
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The total potential $iot( r i ^ 0) = ^ext(?", 6, <fi) + ®ind( r > ®-> 0) for r > i? then becomes 
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On the surface of the sphere, the boundary conditions are $ tot (R, 9, < 



and 



47rfca (R, 9, 4>; w) , where a (R, 6, <fi; to) is the induced surface charge den- 



sity on the spherical shell. Expanding a in terms of spherical harmonics and using linear response theory we find 
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In this equation, 11/, (w) is the SEG polarization function for L an integer and given in terms of the Wigner 3-j symbol 
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where fo(E) is the Fermi-Dirac distribution function and Eg = ft 2 £(i? +1)/ (2m* R 2 ) with £ = 0, 1, 2, • • ■ and m* is the 

(2) 

electron effective mass. The induced potential $\ nd (r; cj) outside the spherical shell may be calculated to be 
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where Sl(oj) = L (ei + £2) + £2 + (2e 2 /i?)Il£(w) is the dielectric function of the SEG. The force on a charge Q at 
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r o = 0"0) 0o> #o) is along the radial direction and can be found using F(r ) = — Q d& ind (r)/d 



yielding 



F(r ) = fcQ 2 ^(L + l) 
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The interaction potential energy £Yi m (ro) may now be calculated from 
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The effective potential is the sum of the image potential and the centrifugal term and is given by [H Q 



^<g(ro,0)=ZC'(r o ) + 



ri 2 x 2 



2M*(r sin6» ) 



(10) 



showing that g5 is not spherically symmetric. In this notation, M* is the mass of the captured charged particle in 

an orbital state with angular momentum quantum number L. In Fig. [5Ja), we calculated V q (ro,0) as a function of 
ro, for chosen L, R. Also, for the background dielectric constant, we chose e\ — 2.4, corresponding to graphite, and 
£2 = 1. The electron effective mass used in calculating the polarization function IIx in Eq. ([B]) was m — 0.25m e where 
m e is the bare electron mass, and the Fermi energy Ep — 0.6 eV. The orbiting particle effective mass M* = m e . 
Figure [2fb) shows how the peak values of the effective potential depend on radius. Of course, the height of the peak is 
linked to the localization of the particle in orbit. In Fig. [31 the ground and three lowest excited state wave functions 
are plotted for the effective potential when L — 2 in Fig. Uta). 

The value of the angular momentum quantum number L as well as the curvature of the surface of these complex 
carbon structures clearly plays a crucial role in shaping the effective potential. Generally, the form for the L-th 
term may be expressed as V^jj (ro) = — Q^l^o 2 ' i+1 ' ) + /Jlt-q -2 , where «l and /3l are due to the image potential and 
centrifugal force, respectively. The coefficient ap is always positive whereas /3l is only negative for L = 0. This 
power-law behavior ensures that no matter what values the two coefficients may have the image term dominates the 
centrifugal term, leading to a local maximum in the effective potential. This is unlike the behavior for a cylindrical 
nanotube where the image term is logaritmic, due to the linear charge distribution, and may be dominated by the 
r^ 2 centrifugal term, leading to a local minimum instead. Consequently, capturing and localizing a charged particle 
by the image potential of spherical conductors and dielectrics is fundamentally different from that for a cylindrical 
nanotube. For the sphere, as shown in Fig. [31 the wave function is more localized around the spherical shell within its 
effective potential, i.e., the wave function is not as extended . Additionally, the confinement of the charged particle is 
close to the spherical surface. 



4 




0.25 
0.20 



0.05 



1.4 



1.6 



R, 





® 




€,=2.4 i 











2.0 



FIG. 2: (Color online) (a) The effective potential V e s between a charged particle and a spherical shell is shown for a number 
of angular momenta L. The radius of the sphere is R — 1 nm and we chose ei = 2.4, £2 = 1. In (b), the height of the peak in 
the effective potential appearing in (a) is plotted as a function of the radius. 




FIG. 3: (Color online) The wave functions for the ground state (n = 1) and first three excited states (n — 2, 3, 4) are plotted 
for the effective potential V c g between a charged particle and a spherical shell when L = 2. We chose e\ = 2.4, ei = 1 and 
R — 1 nm 



The choice for the radius does render some crucial changes in to make a difference in the location and height of 
the peak. However, only the higher-lying localized states are affected. In Fig. [5J we show how the peak height changes 
as L is varied. Additional numerical results corresponding to ei 3> £2 have shown that a spherical metallic shell has a 
reduced potential peak for confining the captured charge. Thus, the spherical metallic shell is not as susceptible for 
particle confinement in highly excited states as the metallic nanotube. [?j-[l0| This indicates that the dimensionality 
plays a non-trivial role in formation of image states and their spatial extension near the surface of the nanometer-size 
graphitic structure. This direct crossover from a one-dimensional to a three-dimensional regime is not determined 
by polarization effects for the structure in the metallic limit since in the limit e% — > 00 in Eq. ([8]). the Hl(oj) term 
makes no contribution. The difference is due entirely to the geometrical shape in the metallic regime where graphitic 
plasmons fail to develop. For finite values of e%, we encounter the regime where excited particles contribute through 
the polarization function ILl(uj) defined in Eq. ©. The behavior of plasmon excitation as a function of angular 
momentum quantum number L for fullcrenes resembles in all respects the long wavelength (q — > 0) limit of carbon 
nanotubes 32] . Furthermore, in the case of the low-frequency n- plasmons in carbon nanotubes, a surface mode 
may develop for large q, due to the difference in the values for the dielectric constants within the graphitic structure 
and the surrounding medium 



Since the polarization function T\l(ui) vanishes identically for L = 0, the attractive part of the effective potential is 
only significantly modified by screening for a fast-rotating electron. This behavior at zero angular momentum differs 
from tubular-shaped image states for single-walled carbon nanotubes which are formed in a potential isolated from 
the tube Q. The large angular momentum image states for spheres may be probed by femtosecond time-resolved 
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FIG. 4: (Color online) Probability density plots for |\E r .t,n(fo)| 2 /'''o when L = 2 and n = 1,2 (upper panel) as well as n — 3, 4 
(lower panel), where n labels the eigenstates, for a spherical shell of radius R — 10 A. The wave function ^L,n(ro) is a solution 
of the one-dimensional Schrodinger equation with effective potential V e g(ro, 9 — it/2) shown in Fig. [2] We chose ei = 2.4 inside 
the ball, whose outline is shown as a thin circle, and £2 = 1 in the surrounding medium. 

photoemission 13]. Our formalism shows that considering photoionization from various levels of Cgo, the Coulomb 
interaction between an external charge and its image is screened by the statically stretched SEG through the dielectric 
function £l(w = 0). The polarization of the medium II ^ which is driven by the electrostatic interaction is generated 
by particle-hole transitions across the Fermi surface. The polarization also determines the Ruderman-Kittel-Kasuya- 
Yosida (RKKY) interaction energy between two magnetic impurities as well as the induced spin density due to a 
magnetic impurity. 

Increasing radius, the position of the peak moves closer to the sphere as i£ _1 /( 2L ). In the absolute units, the 
position of the peak depends as r^' eak ' > ^ R 1 ^ 1 ^ 2 ^. The typical distances from the surface are between 1.3 R and 
1.5 R. Figure [5] (b) demonstrates how the potential peaks (corresponding to the local maximum for the V e s) depend 
on the radius of the buckyball for various angular momentum quantum number L. Clearly, we see that the potential 
peak decreases with increased radius leading us to conclude that confinement is strongest for smaller buckyballs and 
particles with large angular momentum. For the nanotube, increasing L leads to a reduced local minimum in the 
effective potential and the ability to localize the charge Q. Approximately, the curves may be fitted analytically to 
^ l/R for all considered values of L. However, we found that a better fit for L > 5 would be of the form C1/R + C2/R 2 
where c\,C2 are constants. 

For increased e\, i.e., the metallic limit with t\ 3> £2, we have ^ eiL, so that the coefficient 
((2L + l)[e2/ej,] — 1) — > — 1. In the case of dielectric constant t\ ~ 2.4 for the buckyball, the above mentioned 
coefficients lie within the range from —0.4 to —0.9 and decreases with increasing L. Consequently, for the transition 
to the metallic limit, these coefficients are more affected for states with large L. So, we may conclude that t\ has 
little effect on the position and height of the peak in the effective potential. In fact, for the metallic case, the peak 
is observed to be slightly further away from the center (very little difference ~ 1.37 R compared to ~ 1.31 R for 
R = 1 nm). The height of the peak is only slightly decreased in the metallic case (0.24 compared to 0.27 in the case 
of fullerenes). These numbers are provided for fixed L and the unit of energy is the same as that in Fig. [5J 
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Regarding the wave functions and density plots, Figs. 02 and [4] demonstrate the wave function of a bounded 
electron trapped between the infinite hard wall of the sphere and the potential peak. First, we note that we obtained 
qualitatively similar behavior for different values of L, so the electron states corresponding to the potentials with 
different angular momenta are almost the same. We clearly see that the electron wave functions are not exactly 
localized in the "potential well" due to the asymmetry of the boundary conditions, i.e., infinitely high wall on the left 
and the effective potential profile on the right-hand side. The wave functions corresponding to L — are extremely 
delocalized due to the relatively shallow potential. The fact that the effective potential is not spherically symmetric 
means that for arbitrary angle 9, we must solve a three-dimensional Schrodinger equation. However, for trajectories 
parallel to the x — y plane when the angle 8 is a constant of motion, the problem reduces to a quasi-one-dimensional 
Schrodinger equation involving the radial coordinate. In our calculations, we set 6 — ir/2 so that the captured charge 
is moving in the equatorial plane. For this, the centrifugal term is weakest compared to the image potential, but 
still affords us the opportunity to see its effect on localization. The density plots in Fig. @] show that the innermost 
ring is substantially brighter than the outer rings. This is a consequence of the presence of the r^ 2 factor in the 
electron probability function. In contrast, the corresponding plots for the nanotube [Tol ] do not have the innermost 
ring so much brighter than the outer rings because of the fact that the density function in that case depends on 
inverse distance of the charge from the center of the cylinder instead. This is another unusual specific feature of the 
considered geometry and indicates that the captured charge is more strongly localized for the spherical shell closest 
to the surface for the sphere than the cylindrical nanotube. The lowest bound states for Fig. [3] are in the range from 
— 10 to about— 100 meV, with first few excited states lying very close to the ground state energy. Of course, the bound 
state energies may be adjusted by varying the radius R since our calculations have shown that the peak potential 
decreases as a power function with increasing radius. This property allows considerable manipulation of a captured 
electron and its release to a source of holes for recombination and release of a single photon whose frequency and 
polarization are linked to the electron. This single-photon source may have variable frequency with a broad range of 
applications in quantum information where the message is encoded in the number of photons transmitted from node 
to node in an all-optical network. Gate operations are performed by the nodes based on quantum interference effects 
between photons which cannot be identified as being different. The low frequency photons could be in the infrared 
which is most useful range for telecommunications. Another, more general, practical application and technological 
use of such unique quantum states would be to quantum optical metrology of high-accuracy and absolute optical 
measurements. 
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